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Abstract 

A many-body coherent potential approximation (CPA) previously developed 
for the double exchange (DE) model is extended to include coupling to local 
quantum phonons. The Holstein-DE model studied (equal to the Holstein 
model for zero Hund coupling) is considered to be a simple model for the 
colossal magnetoresistance manganites. We concentrate on effects due to the 
quantisation of the phonons, such as the formation of polaron subbands. The 
electronic spectrum and resistivity are investigated for a range of temperature 
and electron-phonon coupling strengths. Good agreement with experiment is 
found for the Curie temperature and resistivity with intermediate electron- 
phonon coupling strength, but phonon quantisation is found not to have a 
significant effect in this coupling regime. 
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I. INTRODUCTION 



In this paper we study the Holstein-DE (double exchange) model 



where i and j are site indices, c| CT and b\ (c icr and 5, ) create (annihilate) an electron of spin 
a and a phonon respectively, Sj is a local spin, <jj = (1/2) cL^o-o-'Cio-' is the electron 
spin operator {a aa > being the Pauli matrix), L? = S- + of is the ^-component of angular 
momentum, n i = J2a Au c iu an d n \ = ■ The parameter Uj is the hopping integral, J > 
is the Hund coupling, h is the Zeeman energy, g is the electron-phonon coupling strength and 
uj is the Einstein phonon energy, if is a model for the colossal magnetoresistance (CMR) 
manganite compounds with the double degeneracy of the conduction band neglected and 
a simplified form assumed for the electron-phonon coupling and phonon dispersion. The 
electron-phonon coupling in Eq. (Jl|) is of the breathing-mode form, i.e. —g' J2i n i%i in the 
classical limit (where Xi is the phonon displacement), but we regard it as an effective Jahn- 
Teller coupling. 

Hamiltonian (|l]) was first studied by Roder et afl, who treated the Hund coupling using a 
mean-field approximation and the electron-phonon coupling using a variational Lang-Firsov 
approximation. The same authors later used a similar method to study a more realistic model 
for CMR systems!. In this paper we treat both the Hund and the electron-phonon coupling 
using an extension of a many-body coherent potential approximation (CPA) previously 
derived for the DE model!'!. The CPA treats the Hund coupling better than mean-field 
theory and has the advantage over Lang-Firsov variational methods that the whole of the 
electronic spectrum can be studied, not just the coherent polaron band near the Fermi 
energy. In the limit of classical spins and phonons Millis et al used dynamical mean-field 
theory (DMFT) to study another more realistic model for CMR materialsi'i. Here however 
we concentrate on the effects of quantisation of the phonons. 

Our approach has many similarities with studies of the Holstein model using DMFT, 
which have been carried out for the classical phononB and empty-band! limits in which 
the model is a one-electron problem. Indeed the standard dynamical CPA is equivalent to 
DMFT for one-electro n pro blems such as the binary alloyi, the DE and Holstein models in 
the empty-band limiti'00, and the DE model with classical local spinsS. DMFT should be 
regarded as the correct extension of the CPA to many-body problems^. For the current 
many-body problem we regard our CPA as an approximate solution of DMFT, or as an 
extrapolation from the one-electron case. The CPA has the advantage of relative analytic 
simplicity, but does not treat the many-body dynamics as well as DMFT, retaining too much 
one-electron character. 

The many-body CPA derived for the finite S DE model in Ref. § and Ref. |] was based on 
Hubbard's scattering correction approximation for the Hubbard modelS Hubbard's approx- 
imation was derived by decoupling Green function equations of motion (EOM) according to 
an alloy analogy in which electrons of one spin are frozen while the propagation of those of 
the opposite spin is considered (within the CPA). Although more modern formulations of 
the CPA exist Hubbard's EOM approach was found to be particularly suitable for extension 
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to the DE model, where the possibility of electrons exchanging (spin) angular momentum 
with local spins complicates the problem. The resulting many-body CPA was exact in the 
atomic limit ty — > and recovered the one-electron CPA/DMFT in the empty-band and 
classical spin (S — > oo) limits. 

In Sec. y we solve the atomic limit of Hamiltonian (J]), and in Sec. |T| we extend our 
many-body CPA to the Holstein-DE model. The properties of the CPA solution are discussed 
in Sec. [TVL and the special case ( J = h — 0) of the Holstein model is considered in Sec. |V]. 
We give a summary in Sec. [VI]. 



II. THE ATOMIC LIMIT 



In the atomic limit — > H is exactly solvable using the canonical transformation 
H ^ H = e u He u where v = Q-nib* — b) (Ref. (we drop site indices). In the presence of 
electron-phonon coupling the phonon potential is of the displaced harmonic oscillator form, 
and the effect of the canonical transformation is to shift the operators to take account of this: 
b i— > b + and c a \— > Xc a where X = exp(g/u(tf — b)). This transformation decouples the 
Hamiltonian H = Hb + Hf into a bosonic component = um b and a fermionic component 
Hf = —JS ■ a — hL z — g 2 /un 2 where g 2 /w =: \u is the binding energy of a polaron. 

The one-electron Green function G a {t) = —i6{t){{c (T {t),c] J }) can be separated into 
fermionic and bosonic traces using the invariance of the trace under cyclic permutations 
and e vX e v = 1, 



G a {t) = -i0(t) 



Ti f {e-? H f} Tr 6 {e-^} 
Tr f {e- pH fcle m f t c a e^f t }Tr b {e- pH >Xie iHbt Xe- iH '> t } 



Tr f {e-P H f} 



Tr b {e-^} 



(2) 



We evaluate the bosonic traces directly and the fermionic traces using the equation of motion 
(EOM) method, and in the energy representation G a (e) = /f^ dt e let G a (t) obtain 



GJe) 



E 



L (2A^6(a;)(6(a;) + l) 



(25 + l)exp(A(26H + l)) Q=± 
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3 r/3u)/2 W 0+ , p -rf3co/2 jxrO- 

v ' rvrr > ° ' ' n n 



e - aJS/2 + ha/2 + Xcu(l + 25 a+ ) + ur 

r/3u;/2 jxrl+ , -r/3u>/2 vxrl- 

^ r y rvrr I r r rv/T 



(3) 



e + a J/2(5 + 1) + ha/2 + Aa;(l + 25 a+ ) + cur 
Here the weight factors 

W% = (S + 1) (n%nl) - aa (S z n a _ a nl) + (1 - 5 a7 ) (S^a + ^) (4a) 
W% = S (n a _ a nl) + aa (S z n\n<) - (1 - 5 ai ) (S^a**) , (4b) 

a, 7 = ±, 5 Q+ = 1 for a = + and for a = 



b(uj) = l/(exp(/3o;) — 1) is the Bose function, and we define n+ = n a , n a 



l r is the modified Bessel function, 

1 — n CT , and 
S-°a +a = S-a+ for a =| and S+a" for a =|. 

In the (J = h = 0) case of the Holstein model Eq. ([$]) reduces to the formula 
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G n 



E 



I r (2A A /6(w)(6(w) + l) 
exp(A(26(cj) + 1)) 



f^ 2 (nZ a n a ) + e- rl3w / 2 (nZ a n-) 
e + Xu + tor 



e + 3\lu 



lot 



(5) 



but we are mostly interested in the strong Hund-coupling limit, so we shift the energy to 
have the zero near the Fermi level, e i— > e — JS/2, and let J — ► oo. In this limit we find 



I r (2\y/b(u)(b(u) + 1)) e ^/2 w + + e -r/3^/2 jp- 
(25 + 1) exp(A(26(u;) + 1)) e + ha/2 + tor + Xcu 

where the weight factors 

W+ = ({3 + 1 + aS z )n a + S-°a+°) = (2S + l)(^ + o (a z )^j 



w: 



((1 - n)(S + 1 + aS z )) =(S + l)(l-n) + a (S z ) - 2Sa (a z ) . 



(6) 



(7a) 
(7b) 



The paramagnetic state spectrum of Eq. is plotted in Fig. [L] for the classical spin limit 
S — > oo at quarter-filling n = 0.5. The spectrum consists of delta-function peaks separated 
in energy by u, and for clarity we include the peaks' envelope curve in Fig. [TJ. Note that the 
symmetry of the spectrum about zero energy is due to choice of filling n = 0.5; in general 
the lower and upper (zero temperature) 'bands' have weights n and 1 — n respectively. By 
counting weights it may be seen that for any n the zero temperature chemical potential 
fi(T = 0) = lies in the zero energy peak in the middle of the pseudogap. 



III. THE CPA GREEN FUNCTION 

We now derive a many-body CPA for the one-electron Green function G l J(e) of the 
full Hamiltonian (p. As discussed in the introduction we proceed by decoupling equa- 
tions of motion (EOM), adapting decoupling approximations previously used for the 
DE modelS. Recall that with the fermionic definition of Green functions, ((A\C)) e = 
—i Jo°°dt exp(iet){{A(t),C}), the EOM is 

e((A;C)) e = ({A,C}) + (([A,H];C)) e . (8) 

As in Ref. f| and in the original version of this method due to Hubbardlil we split the Green 
function into a low-energy component G l J~(e) = ((n~c ia ; c] CT }) e and a high-energy component 
G»+(e) = ({nfc ia ; c} a )) e , i.e. G$(e) = £ a =± G^ a (e). To close the system of EOM we in fact 
need to introduce the Green functions 

Sr(r,e) = ((r i (r)n?c ta ,cl)) e (9a) 
7T (f, e) = ((riCrX-STV,,; c)„)) t (9b) 

where r = (p,(p,9) is a parameter and the operator Ti(r) = exp(pSf) exp(0fe|) exp(9bi). 
Note that T is a generating function for the operators S z , and b, so that d n / 'd(f> n Y 'j(f) = 
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(b\) n Ti(r) for instance. This is convenient in allowing us to close the system of EOM with 
a minimal number of equations. 

When writing the EOM we use the convenient commutation identities [e* 6t , b] = 
— <j)e^ b and [e eb , tf] = 8e eb , and the Feynman operator disentanglement relation e A+B = 
e A e B e -i/2[A,B]^ which hoMs if ^ 5 j ) A j = ^ B} 1 B}= 0. We also work for a =]; the a =| 

equations can be obtained using the symmetry of H. We introduce the operator 



K c 



( 9 d n . ,„ 



(10) 



and obtain the (exact) EOM 



h J d T ^ ni 
" + 2 + 2d-p +K{ 



S\ 3a {r 



Cit: c 



k 



for Sl 3a (r, e) and 



ft J ( r d s 



J 



-pSa 



+ K C 



d d 2 
S{S + l) + a-^- 



nf v H 



dp dp 2 



S l f{r,e 



—a5ij (T(f)S a 



(12) 



for Tl JOI (f, e). Here Ho is the electron hopping term of the Hamiltonian and we have dropped 



site indices in the expectations, assuming a homogeneous state. Equations ( |TTD and ([12]) 
should be compared with their analogues for the DE model case: equations (4) and (5) in 
Ref. |. 



As usual we now neglect the penultimate Green functions in equations (|TT[) and (|12l) (the 
ones containing [n" a , Ho}). This corresponds to making the alloy analogy. We treat the final 
Green functions in these equations using a CPA and treat all other terms exactly. In fact 
we use the approximations 



r(r>?) (E^Gf(e) " Mz)G} j (e) ) + Me)S* a (f,e) (13a) 



^((Urjn-Src 
k 



a (r(f)5-a+) ( E^Gf (e) - J T (e)Gf (e 



+J l ( e + / i )Tf(f,e) (13b) 



where Jo-(e) = e — S CT (e) — Go^e) -1 , S CT (e) and G a (e) being the self-energy and local Green 
function respectively. The function E a (e) = e — J CT (e) is related to the Weiss function of 
DMFT, E a {iuj n ) = Q^iiuJn), iu n being a fermionic Matsubara frequency. Equations (|13a| ) 
and ( |13b| ) are generalisations of Hubbard's scattering correction approximation!^, and lead 
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to the CPA equations in the case of the DE model. We make these particular approxima- 
tions since Eq. ( 13a ) and Eq. (|13b| ) are of the usual CPA form, but do not give a formal 
justification. 

We define E h a {e) = + and \*{e) = 6 tj + E k t ik G k J{e) - J a {e)G»{e), and 

from equations ([0]), (0), ( pL3a| ) and ( |13b| ) obtain 



J_ p p5 a + 
2 C 



2 C 



« Af (e 




ija. 



r, e 



r(r)n? 



a{T(r)S-a + ) 



(14) 



ija. 



We make no further approximations. We use the top row of Eq. fll4j) to eliminate 
obtaining a second-order linear (parabolic) PDE for Sl 3a (r,e). We take i — j and use 
A"(e) = 1 (Ref. |4|). In the strong Hund-coupling limit J — > oo, which we take with the 
energy origin shifted as e i— > e — JS/2, this second-order PDE simplifies to the first-order 
linear PDE 



;i + g)Ef(e) + g£f(£) 
25+1 



£f(e) 



25 + 1 dp 



T(f) 



(S + 1 + S z )nJ + S~a~ 
25+1 



(15) 



Note that in this limit we may assume a = — . We change 
g/uj){cj) + g/oS) and O = l/a;ln(0 + g/u>), in terms of which 



variables to $ 
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K(Q, $) = g ( ^ + $e-" e ) + -- 



90" 



(16) 



In the new (p, 0, <3?) system of variables Eq. (|15D contains derivatives with respect to p and 
only, facilitating its solution. We find (see appendix) for p = 9 = <p = 



Gr(e) 



E 

m,n=0 



m\ n! e A 



x 



e -«M t (6t)«»(6 _ g/oj) n e 9 l ub [{S + 1 + 5 2 )r^ + S~a 
(S+1 + S z )E^e) + (S- S z )E>{(e) + (25 + l)(Acu + (m 



(17) 



where () denotes quantum and statistical averaging and S z in the denominator acts on the 
left. 

In principle the average in Eq. ( |i~T|) should be determined self-consistently, but this is 
difficult to carry out. In previous many-body CPAs, e.g. Hubbard's for the Hubbard model 
and ours for the DE model, it was found that the hopping does not affect the total weight 
in a band near a given atomic limit peak, at least when the bands are separated so that the 
weight associated with a given atomic limit peak is a meaningful quantity. For simplicity 
we therefore assume that all averages take their atomic limit values. Note that owing to the 
degeneracy of the atomic limit states this says nothing about the spin polarisation. This 



6 



assumption means that we cannot take account of the effects of electron hopping on the 
phonon distribution. 

The right-hand side of Eq. (|T7j ) then depends on the half-bandwidth W only through 
E^(e). We change summation variables (m, n) to r = m — n and s = m + n and use local spin 
projection operators P(S Z = m z ) to pull the denominator of Eq. (|I7]) out of the average. 
Since \im w ^ [(S + 1 + m z )E\{e) + (S - m z )Ef(e)) = (2S + l)(e + h/2) we can match our 
averages (summed over s) with atomic limit peak weights to see that 



^ It (2Xy/b(u)(b(u) + 1) 



exp(A(26(w) + 1)) 

/ e- r ^/ 2 (l - n){S + 1 + S z ) + e r M 2 [{S + 1 + g 2 )n T + g-<r+] \ 

\ (5 + l + ^)^(e) + (5-^)^(e) + (25 + l)(Aa; + r^) /' ( ' 

This should be compared with the atomic limit expression Eq. (^j) , to which Eq. ([18]) reduces 
as W -> 0. 

Now in the case of the empty-band limit of the Holstein mo de!0i one is used to obtaining 
a CPA/DMFT expression for G a in the form of a continued fraction. In Eq. fll8|) we have a 
simpler expression involving a sum over the atomic limit peaks, despite the more complex 
nature of the problem which we are considering (the many-electron case with both Holstein 
and DE interactions). One might suspect that our CPA is cruder than the one-electron 
CPA. In fact our expression for G\ for the Holstein model (given in Sec. [Vj) in the limit 
n — ► is equivalent to making the approximation 

E h a {e + ru) w E*{e) + ru (19) 

in the one-electron CPA expression. We will mainly consider the case of an elliptic bare 
density of states (DOS), D(e) = 2/(nW 2 )VW 2 - e 2 , for which it be shown that E(e) = 
e — W 2 G(e)/4. For the elliptic DOS approximation Eq. (|l^) is thus equivalent to neglecting 
energy shifts in the Green function on the right-hand side of the CPA equation. Since we 
do not recover the one-electron CPA/DMFT as n — > 0, unlike in the case of the bare DE 
model@, our CPA for the Holstein-DE model is not as good as our CPA for the DE model. 
We choose however to accept the increased crudeness of the approximation in return for 
the greatly increased simplicity; a CPA which correctly reduced to the one-electron CPA as 
n — » would probably be analytically intractable in the many-body case. 



A. Calculation of Curie temperature 

Using mean-field arguments Millis et a0 claimed that the bare DE model predicts a 
Curie temperature T c at least an order of magnitude too large. However, subsequent more 
reliable treatments of the DE model taking into account quantum fluctuations^^ showed 
that the DE model's Tc is in fact in reasonable agreement with experiment. Now as discussed 
by Roder et afl phonon coupling suppresses Tc- We therefore calculate Tc to see if a phonon 
coupling strength g exists that gives a much larger resistivity (than the g = case) without 
making Tc unphysically small. 
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For simplicity we work in the classical limit S = oo, in which E^(e) = E a (e) since 
h ~ 1/S, and specialise to the case of an elliptic bare DOS, where as mentioned above 
E a (e) = e — W 2 G a (e)/4: is just a function of G a . We set h = in Eq. ( Jl8| ) and expand G a 
about the paramagnetic state to first order (in 5(<r z ) or 5(S Z )), obtaining 



SG^e) = (^2w r (e)e- r ^ /2 ^j 5 (S z ) + 4 (l>r(e) sinh(r/3^/2)^ 5 (a z ) 



(20) 



where 



L (2Xy/b(u)(l + &H) ) (g(e) + ur)- 1 
Wr(e) ~ 2exp(A(l + 26(cu))) 1 + (W 2 /12)G'(e)/E'(e) ' ( } 

A'(e) = dA(e)/de and we drop spin indices on paramagnetic state quantities. From the 
spectral theorem (a z ) = T[G^ — GJ/2 where 

1[A\ = f°° -f(e - n)Im [A(e - z0)] , (22) 

J-oo 7T 

f{e — fi) being the Fermi function. Applying X to Eq. ( |20"D and rearranging leads to 

' (a)= l-4E,XK]sinh(r^/2) 5(5) - (23) 

In Ref. |] we showed (in the DE model case) that the CPA for electronic Green functions 
does not give a good estimate of local spin expectations. Fortunately for S = oo we can 
use DMFT to obtain an expression for S(S Z ) in terms of SG^ which is exact in the infinite- 
dimensional limit. After integrating out the bosonic degrees of freedom the DMFT effective 
action can be written in the Matsubara representation as 

o__ V /i t \( E ^n) + (J/2)S Z (J/2)S- \(CnA_„ h qZ 

n ^ T 1 ' V (J/2)S + E^) - {J/2)S Z ){c nl ) ^ 

+ [ dr [ &r'ni{T)U{T-T')ni{T') (24) 



JO 



where we work for S = oo, h and J finite, and the n subscripts refer to fermionic Matsubara 
frequencies. The last term in Eq. (^4|) is an attractive Hubbard-like term, with the Fourier 
transform of the interaction given by U{iuj n ) = — {l/2)g 2 /{ijj 2 + u 2 ). It is retarded in 
imaginary time and originates from the phonon couplingp. We expand S about the h = 
paramagnetic state with action S and partition function Z , S = S + 5S + ■ ■ ■ where 

5 S = -p h S z - £ (5£ T (^ n )4 T c nT + SE^iu^c^) , (25) 

n 

and in terms of 5S we have 

S (S z ) = ~yJ d " S (n dc Ld Cr J S z e- §0 5S (26) 
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where / d 2 S is the integral over the surface of the unit sphere. 

Now SEaiiujn) = —W 2 5G a {iuj n ) /4 for an elliptic DOS, and we use the relation 
J2 n di^n) — Z[g], which holds for functions g analytic off the real axis, to write Eq. 

(f26|) as 



{ h W 2 
6{SZ)=(3 }3~~ 1 



'dS{r,e) 



5Gt 



r=Q 



(27) 



Note that the effects of phonon coupling enter only implicitly via Green functions. This is 



expected as the electron-phonon coupling is spin-symmetric. Setting h = in Eq. (|27|) and 
using (|20|), (p3|) and dS(f, e)/dp\?=o = G(e)/3 we obtain the Curie temperature equation 
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E r Z[Gw r }e- r W 2 + 4E„I[Gw r ]I[w 8 ] smh(/3u(r - s)/2) 
1 -4£ r Z[uv] sinh(r/^/2) 



(28) 



upon dividing by S(S Z ). We discuss the value of Tq in the next section. 



IV. RESULTS 

We now discuss numerical results obtained using our CPA, for simplicity using the elliptic 
bare DOS and working at J = S = oo and n = 0.5. In the spin-saturated state the minority- 
spin weight at low energy is of order 1/5*, so the classical limit S —>■ oo is convenient as 
we do not need band-shifts, which are difficult to obtain within the CPA, for consistency 
of the saturated ferromagnetic state. Quarter-filling n = 0.5 is used because owing to the 
symmetry of the spectrum about zero energy the chemical potential fi(n = 0.5) = for 
all T. For a homogeneous state the doping has a qualitative effect only as n — > or 1, 
and the results for a more physical value n ~ 0.7 are similar in form to those at n = 0.5. 
Note however that the quantitative predictions of the model are very sensitive to the model 
parameters, especially the electron-phonon coupling strength but also the doping, so if the 
physical doping value is used the model parameters must be adjusted to retain quantitative 
agreement with experiment. We take u/ks — 0. 05^/^6 ~ 600K for W ~ leV. Zhao et al 
report that u/Jcb ~ 100K for Lai_ x Ca x Mn03 (Ref. |17|) so this may be a bit large. 

The Curie temperature Tq obtained from equations (p8|) and (|18|) is plotted against 
electron-phonon coupling strength g in Fig. |2|. It will be found later that g ~ 0.161^ gives 
reasonable values for the resistivity. For values of g in this range Tq is only suppressed by 
about a factor ~ 2, so for W ~leV is still compatible with experiment. 

The effects of phonon coupling on the (forced) h = T = paramagnetic state DOS are 
shown in Fig. |3|. At g = we obtain the usual elliptic bandi. As the coupling g is increased 
the DOS broadens, small subbands are split off from the band-edges, and a pseudogap 
appears near the Fermi energy. At a critical value g c the DOS splits near zero energy, 
leaving a small polaron band in the gap with low weight but very large mass. Increasing 
g further causes more bands to be formed in the gap, with weights equal to the relevant 
atomic limit peak weights. The effect of increasing the temperature T on the DOS in the 
pseudogap is shown in Fig. |] for g = 0.18W > g c . With increasing T the DOS at the Fermi 
surface increases rapidly and the polaron bands are smeared out. 
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For g > g c the majority of electrons are in fully occupied bands, and the itinerant 
electrons lie in a polaron band of very small weight near zero energy. At T = this band 
is equivalent to the one obtained in the standard strong-coupling theory of the Holstein 
model0, where one averages the phonons out of the Hamiltonian H considering only diagonal 
electron hopping processes in which the number of phonons in each state is conserved. In our 
approximation this polaron band is damped even at T = (i.e. ImS(/x) 7^0), but it would 
be coherent (barring damping coming from the disordered local spins, i.e. in the saturated 
ferromagnetic state) in an approximation which took better account of the dynamics. In the 
usual strong-coupling treatment, which treats only the coherent polaron band, it is found 
that the DOS D(fi) at the Fermi surface decreases with increasing T. This is not inconsistent 
with our finding that D(/j,) increases with T since our DOS includes all the spectral weight, 
both (ideally) coherent and incoherent. For n ^ 0.5 the DOS is no longer symmetric about 
zero energy; the main lower and upper bands into which the DOS is split for g > g c have 
approximate weights n and 1 — n respectively. Although these large features of the spectrum 
vary considerably with doping the zero-temperature chemical potential is always confined 
to the polaron band near zero energy, moving from the bottom at n = to the top at n — 1 
(so that we obtain a band insulator in these cases). 

In our CPA we have no reliable means of calculating the probability distribution func- 
tion P(S Z ), so to go below Tq we use the mean-field approximation for the ferromagnetic 
Heisenberg model with classical spins and nearest-neighbour coupling. Since we regard our 
CPA as an approximation to DMFT, which is also exact in the infinite-dimensional limit, 
this simple approximation may not be unreasonable. We obtain the coupling constant for 
the Heisenberg model by matching Curie temperatures. We take this coupling constant to 
be temperature- independent, but in a more systematic mapping onto the Heisenberg model 
one would expect the coupling constant to vary with temperature. Note also that in prin- 
ciple the Heisenberg model's P(S Z ) is of a different form to the DE model'sB. One effect 
of using a mean-field approximation for the magnetisation will be to obtain the mean-field 
magnetisation exponent of 1/2; in three dimensions we expect the magnetisation to increase 
more rapidly below Tq, but note that Schwartz et aJll find that the magnetisation exponent 
in La . 8 Sr . 2 MnO 3 is 0.45±0.05. 

We plot the up- and down-spin DOSs for T = 0.005W/k B < T c and g = 0.16W > 
g c in Fig. |5], also showing the saturated ferromagnetic and paramagnetic state DOSs for 
comparison. The large difference between -Df er ro(/-0 and -D para (/z) mean that for given T we 
expect the paramagnetic state to have a much higher resistivity than a magnetised state. 
Note that there are no separated polaron bands near /i in the up- and down-spin DOSs, 
even at this low temperature k B T = O.lu where quantum effects might be expected to be 
important. The transfer of weight to the up-spin DOS has broadened the polaron bands 
enough to remove the gaps in the DOS, and mixing of the down-spins with the up-spins via 
the Hund coupling suffices to remove the gaps from the down-spin DOS too. It therefore 
appears that the development of magnetisation prevents quantum effects from becoming 
important for this coupling strength, at least as far as the DOS is concerned. 

We calculate the resistivity p using the formula obtained in Ref. ||, plotting it against 
temperature in Fig. |6| for various magnetic fields h. The form of the curve is broadly in 
agreement with experiment!!!, with the resistivity peak (which occurs at Tc) of the correct 
order of magnitude for La .75Ca .25MnO3 (Ref. ^) and we find a large negative magnetore- 
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sistance near the peak. Note that for W ~leV the Curie temperature Tc ~230K and the 
magnetic field B = 0.004W/(<7//b) ~20T (for g ~ 7/2). The main differences between Fig. 
|6] and experiment are the large residual T = resistivity, due to the artificial incoherence 
of the CPA, and the less rapid drop in the resistivity below T c and with h, possibly due to 
the mean-field form used for the magnetisation. The rise in the resistivity below Tc is due 
to the effects of the reduced spin polarisation on the DOS. Below Tc these effects dominate 
over the effects on the DOS of thermal smearing, which is responsible for the fall in the 
resistivity above Te- 
la. Fig. [7] we show the effect on the resistivity of increasing hydrostatic pressure, which we 
model as an increase in the bandwidth. Any change in the other terms of the Hamiltonian 
is neglected (note that the resistivity is proportional to the lattice constant a, so a decrease 
in a will only reinforce the trend observed in Fig. [7]). The strong suppression of the peak 
and the increase in Curie temperature is in agreement with the measurements of Neumeier 
et on Lao.67Ca .33Mn0 3 , where a drop in peak resistivity of a factor of ~ 2 is observed 
when a pressure of 1.62GPa is applied. This change comes mainly from a decease in the 
effective coupling constant g 2 /{uiW). 



V. THE HOLSTEIN MODEL 



We now briefly consider the special case J = h = where Hamiltonian ([[]) reduces to 
the Holstein model. In this case CPA equation (|14|) takes the form 



[E a {e) + K a (9, <P)} ((exp (#J) exp (06*) <c iff ; 4)) f 
and can be solved as in Sec. |TTT| to yield 



exp 



exp (9b) n a _\ , (29) 



GJe) 



E 



I r (2X^/b(u)(b(u) + l) 
exp(A(26(c;) + 1)) 



_ Ejej + \oj + ur 

e T ^/ 2 (n. a n a ) + e~ r ^' 2 {n^ a n-) 
E a (e) + 3Ac<j + Lor 



(30) 



There is now no mixing of up- and down-spins in the problem so our CPA takes the simple 
form (7(e) = G/a,(E(e)) (where Gal is the atomic limit Green function) which one would 
guess for a many-body CPA: Eq. (^) is just the atomic limit result Eq. (^) with e i— > 
E a (e). The first (G° = ~) term of Eq. (|30|) corresponds to polaron bands near energy — \u, 
and the second (G" =+ ) is the bipolaronic term, corresponding to bands near — 3Xu. Our 
approximation's reliance on the atomic limit means that all bipolaron coupling takes place 
on-site. 

In principle we can use the spectral theorem to determine all weights self-consistently 
in terms of the Green functions G", but for low temperature and strong coupling g most 
electrons are bound as bipolarons and we may set (n^n^) ~ (n a ) = n/2. The groundstate of 
the Holstein model is actually believed to be either superconducting (away from half-filling 
and at strong coupling) or a charge density wave (near half-filling and at weak coupling)£3. 
However, determining the weights self-consistently near the homogeneous state we do not 
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find a (second-order) transition to a charge density wave state. This is reminiscent of the 
CPA for the Hubbard model, where no transition to ferromagnetism or antiferromagnetism 
exists. We are also unable to consider superconductivity within our approximation. 

Note that the true CPA/DMFT result in the empty-band limit, first obtained by Sumi0 
using the CPA and rederived by Ciuchi et aB using DMFT, is 



^) = (l-e^)E g(£) _,„ (t) _ B „ (£) (3D 



where A n is the finite continued fraction 

E(e + oj 



2 

Mt) = ~, . ; nQ ^ (32a) 



E(e+2uj) ( "- 2) f 



B(e+nu) 



and B„ is the infinite continued fraction 



{n + l)g 



Bn(e) = — • (32b) 



E(e-iu 



(n+2)g 
~~Ftl oTa (i+3)g 2 ~ 



As mentioned in Sec. [H| our result in the empty-band limit is only equivalent to the true 
one-electron CPA result if the approximation E(e + ruj) ~ E(e) + ruj is made in equations 
O and ( pb| ). 



VI. SUMMARY 



In this paper we have extended our many-body CPA, developed in references || and ^ for 
the DE model, to study the Holstein-DE model, which we regard as a simple model for CMR 
materials. We were interested in effects due to the quantisation of the phonons. Our CPA 
has the advantage over DMFT of being analytically relatively simple, although necessarily 
cruder, and over variational Lang-Firsov approaches of being able to study the whole of the 
spectrum, not just the low-energy coherent polaron band. We solved the Holstein-DE model 
exactly in the atomic limit in which the CPA becomes exact and solved the CPA equations in 
the strong Hund-coupling limit J —>■ oo. Using a DMFT result for the local spin polarisation 
in terms of electronic Green functions we obtained an equation for the Curie temperature Tq. 
For intermediate electron-phonon coupling strength we obtained reasonable agreement with 
experiment for most calculated quantities, including the Curie temperature and resistivity. 
It appears however that for this range of coupling the development of magnetisation below 
Tc prevents the quantisation of the phonons from affecting the DOS near the Fermi surface 
even at low temperatures. 
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APPENDIX: 



We now solve Eq. flT5|) for 5-|- using the method of characteristics. For compactness of 
notation we define the operator U — (5 + 1 + S z )nJ + S~a + /(25 + 1). In the (p, 0, $) 
system of variables Eq. (15 ) takes the form 



(e) - Ef(e) 



25+1 



OSt OS 



+ 



dp oe 



dP s* 0(e)bt o(e,*)6 



e r ' e 

2 



UJ 



+ g§e~ 



25+1 



5 T (Al) 



where 0(0) = e — g/oo and 0(0, $) = $e + g/oo. The characteristic equations are 



dp 
ds 
d5 T 
"ds" 



£f (e) - Ef(e 
25 + 1 



e p5 z e ^e)6t e 0(e, 4 ,)6 f/ \ _ 



CO 



d0 

d7 



-uje 



+ 



1 + S)£f(e) + SEfje) 
25+1 



St. 



The first two are solved immediately as 
P = Po + 



r^(e) - Ef(e) 



25+1 







(A2a) 
(A2b) 

(A3) 



where po is an arbitrary constant and we set the constant in the G equation to zero without 
loss of generality. These solutions are substituted into Eq. ( |A2b| ) , which may then be written 

as 

d 



ds 



X[(l+S)E'l(e)+SE'l(e)}/(2S+l)+ g y^)s-^e- 



St 



3 ( l(l+S)Ef (e)+SE^(e)} /(2S+l)+g 2 /«) a 



x < e 



(po+s[E^-E h l ]/(2S+l))S^ e ~^ e e^ e ^e-^ (b-^^bjj 



(A4) 



We expand exp(e^ s 6^) and exp($e UJS (b — g/co)) in Eq. (|A4j) as series and integrate to find 



5V exp 



;i + S)E* + SE'l g 



25+1 



+ — s - — $e~ 



5 



UJ 



to 



m,ra=0 



mm! 



E £o(^ oS v^ (6t)' 



(25 + i) e ([( 1 + s + sz ) s T+( s " sz ) £ |]/( 2S ' +1 ) +92 / a;+ ( m - m ) a; ) s 
1 + 5 + 5 2 )E T h + (5 - 5 2 )£f + (25 + l)(0 2 /u; + (m - m)u;) 



b _ £. ) e^f/ ) . (A5) 



We then write the characteristics as intersections of the surfaces 5| = Sf (p, O, $) and p = 
p (p, 0, $), and the general solution of Eq. (O) is of the form 5| (p, 0, $) = F(p (p, 0, $)) 
where F is an arbitrary function. Rearranging we obtain 



5t = exp ( — $e 



-cjG 



exp 



1 + S)Ef + 5£f 
25 + 1 



+ \u 



F p 



t5 e 

25+1 , 



+ E 



,(m-n)ii)9 



m,n=0 



m!n! 



(25 + l)e pS2 e-^ 6t (6 t ) m (6 - g/u) n e 9 f" h U 
(5 + 1 + 5*)^ + (5 - S Z )E\ + (25 + 1)(Acj + (m - n)u;) 



• (A6) 



13 



Now from definition Eq. (|9a|) of it may be seen that is of the form 

\ S oo 



Tph rph 



P+ oc ■ i e » e » $ = E E « m n(p,$)exp 
V + i / m= _ 5n= _ 00 



- E\ \ 

m — '— + mo 6 

25+1 / 



(A7) 



The final term of Eq. ( |A6| ) is compatible with this form but the term proportional to F is 
not, so we must have F = 0. Finally, in our original (p, 8, <fi) system of variables 



Si{p,0,<j>) = exp 



6- 



10 



E 

m,n=0 



m! n! 

m 



{2S + l)e^"e-»/ w6t (V) m (&-#/cj) n e^ 6 
1 + 5 + S 2 )^ + (5 - S^Ef + (25 + l)(# 2 /cu + (m - n)u;) 



17 . (A8) 
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FIGURES 



0.03 



(0 

2 0.02 

i 

CO 
CD 
Q. 



- 0.01 

E 
o 

3 



0.00 



k B T/W=0 
k a T/W=1 



Weight n=0.5 



Weight (1-n)=0.5 




-1 1 
Energy (W) 



FIG. 1. Peak weights of the atomic limit spectrum at low and high temperature. Plot is for 
the (forced) paramagnetic state with S = J = oo, h = 0, n = 0.5, tu/W = 0.05 and g/W = 0.18, 
W being an energy parameter. 
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FIG. 2. Suppression of the Curie temperature Tq with increasing electron-phonon coupling 
g/W. Plot is for S = J = oo, h = 0, n = 0.5 and uj/W = 0.05. 



16 



0.5 



g/W=0 
g/W= 0.1 
g/W= 0.15 
g/W= 0.18 




-1.0 -0.5 0.0 0.5 
Energy (W) 



1.0 1.5 



FIG. 3. Development of the zero-energy pseudogap and subbands in the spectrum with increas- 
ing electron-phonon coupling g/W. Plot is for the (forced) zero-temperature paramagnetic state 
with S = J = do, h = 0, n = 0.5 and uj/W = 0.05. 
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FIG. 4. Evolution of the subbands in the pseudogap (see Fig. ||) with temperature. Plot is for 
the (forced) paramagnetic state with S = J = oo, h = 0, n = 0.5, uj/W = 0.05 and g/W = 0.18. 
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FIG. 5. Comparison of low-temperature spectra for different magnetisations: saturated ferro- 
magnetism and (forced) paramagnetism for T = versus calculated up- and down-spin spectra for 
k B T = 0.005W < k B T c where {S z ) = 0.915. Plot is for S = J = oo, h = 0, n = 0.5, uj/W = 0.05 
and g/W = 0.16. 
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FIG. 6. Resistivity versus temperature for S = J = oo, n = 0.5, uj/W = 0.05 g/W = 0.16 and 
various h. The lattice constant a = 5fA. 
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FIG. 7. Effect of pressure (increasing half-bandwidth W) on the resistivity. Plot is for 
S = J = oo, n = 0.5, lo/W = 0.05 and g/W = 0.16. The lattice constant a = 5fA. 



19 



